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Abstract 

The central limit for the products of identically-distributed positive free random 
variables are studied by evaluating all the moments of the limit distribution. If the 
logarithm of the free random variables are zero mean with unity variance, the central 
limit is found to confine between exp(±co) with cq 2.08. This central limit distribution 
is not the same but very close to the log-semicircle distribution. 

Let Xi denote identically-distributed positive free random variables, if \ogXi ex- 
ist and are zero mean, the central limit distribution of y = lim„_>oo IliLi ^/^^ is 
bounded by | logy| < cq. All the moments of Y are found and high-order moments are 
asymptotically the same as those of log-semicircle distribution. If the variance of log Xi 
is small, the central limit is very close to log-semicircle distribution, showing similarity 
to the relationship between normal and log-normal distribution in commuting random 
variables. If the variance of logA^ is large, the central limit of logF is very close to 
uniform distribution. 

Keywords: Free central limit theorem, products of free random variables, log-semicircle 
distribution 
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1 Introduction 



Free random variables are random infinite-dimensional linear operators that are equiy - 
alently very large rand om matrices. Free probability theory Nica and Speichei ( 20061 ): 
Voiculescu et al.l (jl992l ) is used to study free random variables. The t h eory for free ran- 
dom variables are applicable to study random matrices Speichei ( 20031 ) ( Voiculescu et al. . 
1992I . §4). The statistical properties of free random variables are equivalently that of the 
eigenvalues of large random matrices. 

In free probability theory, the central limi t theorem on the sum of independent free 
random variables gives semicircle distribution Voiculescu (1986, 199ll ) ( Voiculescu et al. . 
1992I . §3.5). Semicircle distribution serves the same function as the Gaussian or normal 



distribution for the sum of independent commuting random variables. Mathematically, if 
Xi , X2 , . • • , Xn are identically distributed zero mean free random variables with variance of 
(ii/2)^, the free summation or additive free convolution of 



Xi ffl X2 



Xr, 



n 



has the semicircle distribution of 



p{t) 



1 





- t2 



\t\ < R 
otherwise 



(1) 



(2) 



1 



where R is the radius of the distribution. 

For positive conventional commuting random variables, the product of positive random 
variables has a central limit as the log-normal distribution because of xiX2 = exp(logxi + 
log X2) for positive random variables of xi and X2- As the central limit for the logarithm of a 
random variable is normal distribution, the product of random variables has its central limit 
as l og-normal distr i butio n. Log-normal distribution is applic able to many physical phenom- 
ena iLimuert et aD (|2nnih . the shadowing of wir eles s chan nel (jRappaportl . l2002i . §4.9.2), and 
some nonlinear noise in fiber communi cations iHol (l200dl . L o g- normal distrib ution is also 
very close to the power law or Zipf law iMitzenmached (j2004l ) ; iNewmanI ()2005l ) . 

For positive free random variables, the relationship of Xi M X2 cannot guarantee the 
same as exp(logXiffllogX2). Alternatively, the relationship of exp(XifflX2) is not necessary 
equal to the free multiplication of exp(Xi) Klexp(X2). The log-semicircle distribution is not 
necessary the central limit for the product of positive free random variables. 

Here, the central limit distribution is derived for the products of positive free random 
variables. The moments of the central limit distribution are not the same as that for log- 
semicircle distribution but numerically very close for the high order moments. The central 
limit for the sum of zero mean unity variance free random variables has a radius of 2 as 
from ([2]). Let Xi denote identically-distributed positive free random variables, if \ogXi are 
zero mean with unity variance, the central limit of products of X^^ with n approaching 
infinity has a radius of approximately cq = 2.0805. Although very close, the central limit of 
the product of free random variables is not log-semicircle distribution. Both the difference 
and similarity are presented in later parts of this paper. 



Th ere was a long history to study the products of random matrices iFurstenberg and Kesten 
to represent a physical system with the cascade of a chain of basic random ele- 
ments with identical statistical properties. The results here are bas ed on the multipli- 
cation of free random varia b les IVoiculescul (119871). Similar works of iBercovici and Fata 



Tucci 



(|2nOn|') : IBercovici and Wand (120081 ): IChistvakov and G^tzel (|2008l ): iKarginI (l2007bh 

(|2010l ) were mostly on the behavior of the products of free random variables but none 

of them considered the central limit distribution. The behavior of the products of free 

I ~ 1 J T\ I — ' ' ' 

random variables was studied in IBercovici and Fatal (|2000l ): IBercovici and Wang (2008); 

Chistvakov and Gotze ( 20081 ) as infinitely divisible free random variables. In Kargini ((2007b ) . 



the norm was consider ed for the c ase with unity mean finite support Xi. The geomet- 
ric mean was found in iTuccil (|20in|). similar to the la r ge nu r nber theory for comrnutativ e 
random variables. The theory in Bercovici and Fatal (l200nl ): IBercovici and Wand (|2008l ): 
Chistvakov and Gotzd (|20oi ) : iKarginl (j2007bl ) : iTuccil (|20ld ) may applicable to the case that 



log Xi do not exist, i.e., with an atom in zero. The theory here is only a pplicable to the case 
with t he existing of logXj. T he results here are very siin ilar to that of IBercovici and Fata 
(|2nonl ) but cannot derive from IBercovici and Fatal (|200nl ) directly. 

Later parts of this paper are organized as following: Sect. [2] summarizes the main results 
of this paper; Sect. Ogives the proof of the theorems related to the products of free random 
variables; Sect. [H shows numerical computation of the products of free random variables, 
mostly to verify that the central limit distribution is very close to log-semicircle distribution 
if \ogXi have unity variance; Sect. [5] is the conclusion of this paper. 
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2 Main Results 



The goal here is to find out the central limit for the products of free random variables. 
Although we are not able to write down the limit distribution explicitly, all moments of the 
central limit can be derived. 

Theorem 1. Let Xi denote identically-distributed positive free random variables, if log Xi 
are zero mean and have unity variance, the limit of 

1 1 1 

Y = lim Xf^ M Xf^ Xf^ (3) 

n— s>oo 

has the following properties: 

(a) The k-th moment ofY is 

k /2 

^(yfc) = £_4i)i(-fe). (4) 

(b) The asymptotic expression of (piY^) for large k is 



cok 

with Co = 2.0805 and oq = 1.0034. 
(c) Y is bounded by the relationship of \ logY\ < cq. 



2ao^ (5) 



In Theorem \T\ L^a^ {z) is the Laguerre polynomial Koornwinder et al. I (I2OIOI) (ISzegq . 



1975 . §5.1) and In{x) is the modified Bessel function of the first kind (jOlver and Maximonl . 
2010l . §10.25). 

If log Xi do not have unity variance, we have the following corollary: 
Corollary 1. If the variance of log Xi is a^, the k-th moment of the central limit Y is 

V{Y') = -^41,{-alk), (6) 

and I logy I < Co with 



CO = a,\ll + ^ + 21og ^ + y 1 + ^ . (7) 



The asymptotic of bound cq is cq ~ 2ay for small and cq ~ (t^/2 for large (T„. 

This study is always limited to the case with zero mean log Xi . If log Xi have positive 
mean, Y ([3]) grows to infinity. If log Xi have negative mean, Y ([3]) shrinks to zero. Only the 
cases with zero mean logXj is interested for the study of y ([3]). In practice, the mean of 



log Y may be studied separately as the Lyapunov exponent as in iFurstenberg and Kesten 
(| 19601 ). 



Definition 1. A positive random variable Z is log- semicircle distributed with radius of R 
if log Z has semicircle distribution with radius of R. 
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Prom Theorem [H all moments of the central limit Y is known. In principle, if all mo- 
ments of a probability measure is known, the probability measure can be determined, in 
this case, the H ausdorff moment problems with boun ded interval between exp(ibco) for Y 
( Akhiezer . 19651 . §2.6.4) ( Shohat and Tamarkin . 1943 . §1.2.d). The probability distribution 
can b e uniquely determined becau se the moments are conformed to the Carleman's condi- 



tion ( Shohat and Tamarkin . 19431 . §1.6). However, practical evaluation of the probability 



distribution is difficult. The moments of Theorem [T] can be compared with the moments 
of log-semicircle distribution with radius of cq. When the variance of logXj is small, the 
distribution of Y is very close to log-semicircle distribution. 

The moments of log Y can also be found and used to determine log Y as the Hausdorff 
moment problems with bounded interval between =bco. If logXi have small variance, the 
moments of logY are determined by Catalan number, similar to that of semicircle distribu- 
tion. If log Xi have large variance, the moments of log Y are the same as those of zero mean 
uniform distribution in the interval of ibo"^/2. The distribution of logY is approximately 
semicircle distribution if <C 16. The distribution of logY is approximately uniform 
distribution if cr,^ » 16. 



3 Products of Free Random Variables 

Instead of consider the product of positive identically distributed free random variables Xi, 
we assume that Vi = logXj exist such that ip{v.i) = E{vi} = and ip{vf) = E{vf} = 1, 
i = 1, . . . ,n. The free random variables of exp(t;j), i = 1, . . . ,n, are always positive. Define 

r„ = exp(^)sexp(3=)«...aexpfe), (8) 



and 

Y = lim Yn. (9) 

n— >oo 

Same as ([3]), the free random variable Y is considered here. 



Lemma 1. IVoiculescul (jl987l ) For free random variable A with probability measure fJ.A{t), 
define its moment function as 

and its S-transform as 

Sa{z) = ^XAiz), (11) 

where Xa{z) is the inverse function by composition of the moment function ij) ■ If A and 
B are positive free random variables, 

Samb{z) = Sa{z)Sb{z). (12) 



Using Lemmadlfor Yn ([9]), 5'y„(z) = S^^j^{zY and Sy{z) = lim n-»oo gy^ (z). Although 
Lemma [1] may extend to free random variables that are not positive iRai Rao and Speicher 
(|2003), i'ree random variables exp {vi/^/n) here are always positive. 
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Lemma 2. The S-transform ofY is 



Sy{z) = exp ( -z - ^ ) . (13) 



Proof. We first derive S^y-/^{z). The S-transform of exp {vi/y/n) can be found as a power 
series of l/\/ra, in essence the perturbation analysis of the S-transform. Assume that 
tpviiz) = + X^fcS vl^f )^'^ Vi are zero mean, unity variance free random variables, 
where V'(ff ) = } is the /cth moments of Vi. Note that the moments of tp{v^), A; > 3, is 
not very important to arrive with the results but includes here for completeness. 
The moment function of exp(f j) is 



, /X z z{l + z) ^ ip{v^) 

^ ' k>3 



(14) 



t=o 



where 



t=o 



^i^i I z is the coefficient corresponding to v'(vf ) for the mo- 
ment function ^pe^i{z). The summation of Ylii>i ^'^^^ be obtained by expanding e^""^ in 

The moment function of exp (vi/y/n) as a power series of l/\/ra is 



1-z 



1 + 



1 + z 



,93{z) 



where 



2n(l 

93(2) 
9Akz) 



"I I 3 
6n2 



24n2 



1 + 4z + z2 

z^ + ll^^ + llz + l 



(1 - zY 

The inverse function of (jlSp can also be expressed as a power series of 1/ \/n to 



1 + Z 



n2 



-h3{z) 



hi{z) +0{rr2) 



where 



hz{z) 
hi{z) 



Z'' + Z + 



z-^ + 2z^ + -z + 



The S-transform of e"'/^ is 



S^.jV^{z) = l-^(z+]^ 



24 



^/.3(.) + #/.4(.) + 0(n-t). 



n2 



(15) 



(16) 



(17) 



The S-transform of 1^ ([8]) becomes a power series of l/\/n using the relationship of 

Sy„{z) = [S^.jvrAz)T as 



SySz) = exp 



n 



n 



Lp{vf)h4{z) 



1 1 

2 V 2 



+ 0(n- 



(18) 
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and 



Sy{z) 



lim SySz) 

□ 



exp 



(19) 
□ 



Th e S-transform of (1191) is inf in itely divisible from th e theory of iBercovici and Pata 
(|200d ): IBercovici and Wand (|2008l ): IChistvakov and Gotzd (jioOS). In practice, all terms 
higher than the second-order is not required in (fT5|) , (fT6]) and (fTSt) . The S-transform of (fTSl) 
shows the convergent properties with the increase of n, similar to the analysis of iKargin 
l^'or non-symmetric free random variable Vi with non-zero skewness proportional 
to <^{vf), the convergence is far slower than the symmetric free random variable with zero 
skewness. 

With the S-transform of Y , the moment function of Y can be founded. 

Lemma 3 (Lagrange-Biirmann formula). // the dependence between w and z is given by 
w/(j){w) = z and 0(0) 7^ 0, the inversion series for w = g{z) is given by 



fc>0 



(20) 



This lemma also called Lagrange inversion formula an d was given in ([Whittaker and Watsonl . 

I927I . §7.32). A variation of this lemma was used in iKarginI (|2007al lbl). 



Part a of Theorem\^ The inversion of the moment function iPy{z) is given by 



Xy(z) 



l + z 



exp 



(21) 



Compared with the Lagrange-Biirmann formula of ([20]), (t){w) = (1 -|- w) exp {w + ^2). 
From the definition of moment function in Lemma [H the A;-th moment of Y is given by 



or 



k\ dw^' 



-{l + w 



,k „kw 



w=0 



■ exp 



E 

m=0 



k 



ml \ m + 1 



(22) 



(23) 



by expanding e^'^ and collecting all terms with ^ in (1 -|- w)^e^'^ . 

Using the generalized Laguerre polynomial l!^\z) Koornwinder et al. (l2ninl ) (jSzegol . 



19751 . §5.1), the moment ([23]) can be expressed as 



k 
□ 



4U-k). 



(24) 
□ 



Confluent hyper geometric function iFi{a;h;z) (jSzegol . Il975l . §5.3) can also be used in- 
stead of the Laguerre polynomial. The moment (|24|) becomes 

(25) 



^piY^) = e''/\Fi{-k + 1;2; -k) = i{-k)/k, 
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where M^^^{z) is the Whittaker function (jWhittaker and Watsonl . 119271 . §16.1 



The asymptotic expression for large k for (|24p determines the upper bou nd of Y . The 
Laguerre polynomial in (|24p has different simple asymptotic expressions (jSzeg i Il97,i 

(|2008l l. The expression 



58.22 



that is for L^n\- 



< X ^ n as from iBorwein et al 



( Szegol . Theorem 8.22.4) may extend to negative number to give rough asymptotic 

—x) ^ -y^/n/x/ i (2^/n x), where Ii{x) is the modified Bessel function of the first 
kind ( Olver and Maximon . 20ld . §10.25). Ignore the smooth term, the rough asymptotic 
gives ^(Y^) ~ -^e^^ /yf^k by the asymptotic of Ii{x) ~ /yj2'Kx ( Olver and MaximorJ . 
2010l . §10.40.1). Although for < X <C n, the simple asymptotic gives very close result for 
the case n 



X. 



l + 0{k-% 



Lemma 4 (Asymptotic of Laguerre polynomial). For large k and /3 > 0, 
where 

c = 




(26) 



(27) 



The asymptotic appro ximation (j26p is a special case that can be derived directly from 
([Frenzen and Wond . Il988l . Eq. 4.7) with higher-order te r ms th ere. The s ame asy i nptoti c 
approximation can also be found in (|Koornwinder et al.l . l20ld S18.15.19) iTemmel (ll99fll V 
Using the asymptotic approximation of Whittaker function ( Olde Daalhuis . 2O10l . S13.21.6). 
the asymptotic expression ()26p can also be derived. 



(28) 



Part h of Theorem d Using Lemma [J] for ()24p , part b of Theorem [T] is with /3 = 1 in (|26p , 
we obtain 

-^-2ao^[l + 0(fc-)], 

CqK 



with Co = \/5/2 + 21og (l/2 + ^/5/2) = 2.0805 and qq = c^^ x 5-3/2 = 1.0034. □ □ 

Note that the above evaluation is very close to that using (|Szegol . ll975l . Theorem 8.22.4) 
that gives cq = 2. 

Lemma 5 (Markov inequality). Let X be positive random variable and h as a non-negative 
real function, then 

Pr (hiX) >a)< (29) 
a 

for all a > 0. 



The format of Markov inequality like Lemma[5]is the same as that in (jGrimmett and Stirzakerl . 



1991 §7.3.1). 



Proof of part c of TheoremUl For positive Y and r, and h{Y) = Y^, Markov inequality 
may be written as 



Pr {Y >r) < 



(30) 
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With k — )• oo, using the asymptotic of Ii{z) ~ e^/\/27rz ( Olver and Maximon . 20ld . 
§10.40.1) for dM]), we obtain 



Pr (Y > r) = if r > exp(co) = hm (p(Y 

k—)-oo 



k\l/k 



where cq is the same as that in ()28p . 
With 



1 



hm Xn 

and fohowing the procedure from ()14p to ()3ip . we may obtain 

Pr (y^-*^ > r) = if r > exp(co 



]_ 



(31) 



(32) 



(33) 



The only difference between Y and is that sign of all odd-order terms in ()14p is 
negative. All the odd-order terms do not change the results to (|3ip . 

With both ([3T]l and ([33]) . y is bounded between exp(ibco). Equivalently, |logy| < 
Co. □ □ 



With all the moments of the central limit of Y from (j24p . the probability measure is 
determined as a Hausdorff moment problem. The Hausdorff moment problem has a solution 
if ^iY^) exp(— cqA:) is completely monotonic ( Akhiezerl . ll965l . §2.6.4) . We numerically verify 
that v^(y^) exp(— CqA;) is completely monotonic up to A: = 200. 

Remark 2 (Saddle point approximation for parts b and c of Theorem [T]). The asymptotic 
expression of 128\) may be found using saddle point approximation to sum each term in Ii23\). 
For large k, binomial coefficient is 



k 

m + 1 



y/7rk/2 



exp 



(m + 1 - k/2y 
k/2 



(34) 



that is based of the normal approximation o f binomial distribution by de Moivre- Laplace 
limit theorem I'Grimmett and Stirzaker. 



199M , UU J ). Th e f actorial ofm\ ~ A/2^ m"^e-"" 



from the Stirling formula Diaconis and FreedmaA i lOSd ). Iwhittaker and Watson . 192 
112.33). 

For large k, the moments [2^) can be rewritten as 



irm 



^exp(-/(m)) 



with 



f{m) 



, (m+l-A;/2)2 
iog^Kj * m H — h mlogym) — m. 



kl2 



The saddle point is located at tuq = nk with 



1, 



4exp 2 



k 



-Wq (46^) ^ 0.62 



(35) 
(36) 

(37) 



and f"(mn) = 4:/k + 1/mo, where Wq is the Lambert W function at the main branch 



Corless et al\ 1199 A ). 
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We may obtain 



kJ27r{K + l)k 



:exp(coA;), 



where 



Co = log 2 + K + 2k^ 



2.0807. 



(38) 

(39) 
□ 



The bound of cq in (|39p has almost no difference with that obtain based on Lemma H] 
although the approximation ()34p is valid only when the saddle point rriQ is in the middle 
and close to k/2. 

If the variance of log Xi is not equal to unity, Corollary [1] may be used. 



CorollaryUi With (p{vf) = a^, following the procedure from (fH|l to (fT9|) . we get 

1^ " 



Sy{z) 



exp 



cr„\ z + 



Following the procedure from (|2ip to ()24p . we obtain 



(40) 



k 



LtU-alk). 



(41) 

with Co given by C, 

3 

• (42) 

2(j|(4 + a2)i 

Following the procedure from ()30p to ()33p . we obtain | logy| < co with co given by ()27p or 



The asymptotic of Laguerre polynomial of Lemma H] has /3 = o"^ in 
([271) and 



Co = C7„'\/l + ^+21og 



□ 



(43) 



□ 



The relationship of ao versus is helpful to understand the asymptotic (j4ip as com- 
pared with that with cr^ = 1. For small cr^ ^ 1, ao — ?• 1 and similar to the case of cr^ = 1. 
For large » 1, ao ~ \/2(Tt,/8 proportional to but increase far slower than cr^. Except 
if cr^ ^> 32, we may assume that ao is close to unity. 

The 2A;-th mo ments of semici rcle dis tribution is given by the k-ih. Catalan number 
Ck = 2FPT(fc') as (|Bressoudl . [2OI0I . §26.5) (jNica and Speicherl . [2OO6I . §2) 



2k 



(44) 



Lemma 6 (Moments for log-semicircle distribution). Let Z as positive random variable 
having log-semicircle distribution with radius of R, the k-th moment of Z is 



viz' 



.him 

kR 



(45) 
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Proof. From definition and using ([2]), the k-th moment of Z is 

<-+R 



^(Z'') = —r / e'^V^^ - 2;2dx. (46) 



With X = RsinO, the integral becomes 

2 /-"/^ 



ip{Z'') = - r e''^''''''^cos'^ed9. (47) 

First move the integration i nterva l to to n, the int egration becomes loikR) — l2{kR) by 
using dOlver and Maximonl . [ioiol . §10.32.3). Using (jOlver and Maximoni . lioiol . §10.29.1), 



loikR) - hikR) = Ii{kR)/kR and 

^(Z') = 2^^. (48) 

□ □ 

Comparing (j45p of log-semicircle distribution with the asymptotic of (j28p for 93(y*^), the 
high-order moments of Y are very close to those of log-semicircle distribution when log Xi 
have unity or less variance. The central limit for the product of free random variables Y ([3]) 
has similar high-order moments as those for log-semicircle distribution with R = cqvo. ()45p . 
This similarity requires that cxq of (|28p is close to unity. Even for the case when log Xi have 
large variance, the high-order moments of Z are close to those of Y because oq in (j42p just 
increases with instead of . 

In additional to compare the moments of (|24p to (|45p , the moments of log Y may compare 
with those of a semicircle distribution given by (|44p . Using (|25p for arbitrary cr^, the 
exponential moment generating function of log Y is 

(^(e^iog^-) = ^(y«) = e'^l^'\F^ (1 - s- 2; -als) . (49) 

Using the Kummer transformation (jOlde Daalhuid . lioiol . §13.2.39) with 



e-^/2^Fi (1 - s; 2; -a^s) = e-^^^/^Fi (l + s; 2; a^s) , (50) 

the moment generating function ()49p is an even function with <^{Y^) = (p{Y^^). All the 
odd moments of log Y are equal to zero. For any positive integer k, the 2/c-th moments of 
logy is given by 

(^(log^^- y) = lim e-"^/\Fi (1 - s; 2; -a^s) . (51) 



The expansion of 

oo 

iFi{l-s;2;-a^s) = 1 + 



ds'^'' 

af{s-l){s-2)---{s-k)s'' 



k\{k + l)\ 

l + ^fks' (52) 



k=l 

oo 



gives the coefficient as 



k=l 



k = iJk= 2^ ^iTTTTv ^^^^ 
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where s{n,k) is the Stirhng number of the first kind ( Bressoud . 20ld . §26.8), and \k/2] 
the smallest integer larger than or equal to k/2. The moments becomes 

2*: f Ak-2i 



IS 



For some low-order moments, we have: 



i=0 



{2k 



Corollary 2. If the variance of log Xi is a^, we have: 

4 

^^^ + Y + 80' 



v9(iog^y) 
^(log^y) 



5af, 23a. 
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+ 



fJ, 



12 



448 



(54) 



(55) 
(56) 
(57) 



to cjf . The 



The 2fc-th moment of f{log'^^ Y) is always a linear combination of ci^'^ 
term dominates for small cj^ but the a^^ term dominates for large cr^. 

The coefficient of o"^'^ for (/^(log^'^ y) is the Catalan number C^. The Catalan number is 
given by z = 2A; in ()54p using /2fc that fu rther uses the lowest order term of i = k in ()53p 
with s(fc + 1, /c + 1) = 1 (jBressoudl . I2OI0I . §26.8.1). 

The coefficient of trf for (/^(log^^ y) is given by 



2k 



(2fe)!E 



s(i + l,l) 



1 



=0 



22fc-«(2fc-i)!z!(i + l)! (2fc + l)22fc 



using s(i + 1,1) = {-lyU teressoudl . lioinl . §26.8.14) and ESo(-2)* 

If cj^ is small, the first term of (^(log^'^ Y) dominates, we get (/^(log^ Y) 
pared with ()44p . logy has semicircle distribution with radius of i? = 2(t„. 

j£ ^2 — +1 — 1„„+ J-™™ „-P ,.,('1 — ,2k 



-- 1. 



is very big, the last term of (p{log 

^(iog2*=y) 



' y) dominates, we get 

2k 



1 



2k + 1 



(58) 



Com- 



(59) 



that is the 2/c-th moment of a uniform distribution between ±(7^/2. 

Because moments of (j59p must be far larger than C^o"^^, with ~ 4^/ V irk^ (Bressoud, 

2010 . §26.5.6), the variance o"^ needs to be far larger than 2^ for ()59p to dominate. Con- 
versely, £7^ needs to be far smaller than 2^ for (j59p to be negligible. 



Corollary 3. // i/ie variance of log Xi is al and <C 16, the distribution oflogY with Y 
given by (0) is approximately semicircle distribution with radius of 2a^ . 

Corollary 4. // the variance of log Xi is and ^ 16, the distribution of log Y with Y 
given by ^ is approximately uniform distribution between ±(7^/2. 

Using large matrices for free random variables. Corollary S] is for matrices Xi and Y 
havin g a very l arge co ndition number. 

In iKarginI (l2007bl). Xj has unity mean with ip{Xi) = 1. Using Jensen inequality 
( Grimmett and Stirzaker . 19921 . §5.6.15) for convex function — log(x), c/?(logXj) < logc/?(Xj). 
If Xi has unity mean, the mean of logXj is less or equal to zero and the products of Xi 
may shrink to zero. 
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Figure 1: Simulated central limit distribution of log Y for cr^ = 1 and = 16 and compared 
with semicircle distribution. The x-axis is normalized by cq (03]). 

4 Numerical Comparison 

From the above results, it may be very certain that 



In the limit of n — t- oo, the left hand side of (|6Up is a log-semicircle distribution with radius 
of 2c7^. Although the logarithm of right hand side of ([60]) or \ogY has an upper bounded 
of Co, this upper bound is determined by (jl3]) that is typically larger than 2(T^. Of course, 
the distribution of Y is numerically very close to log-semicircle distribution when cr^ < 16. 

The simulated distribution of log Y is shown in Fig. [1] when the variance of log Xi 
is either cr^ = 1 or 16. Each free random variables of Vi (or logXj) are approximated 
by 512 X 512 positively defined complex random matrices with eigenvalues determined by 
unbounded normal distribution with the variance of either = 1 or 16. The distribution 
of Xj i s log-n ormal distribution that does not have an unique solution from its moments 
Hevdd (|l963l l. The products of Yn is for n = 256 random matrices. The multiplication of 

1 /2 1 /2 

X1X2, for example, is given by X{'"X2X{'\ The probability density of Fig. [U is given by 
the histogram of 51,200 eigenvalues. 

From Fig. (T] the eigenvalue of log Y256 is very close to semicircle distribution if log Xi 
have unity variance of o"^ = 1. The discrepancy to semicircle distribution is very difficult 
to observe. With large variance of cr^ = 16, log 1256 h^is observable difference with respect 
to semicircle distribution but still very close. With a1 = 16, the distribution of log 1256 is 
almost the same as semicircle distribution near =bco. The probability density has the trend 
to be more uniformly distributed and is larger than semicircle distribution in the third 
quarter point of ±0.75co but smaller than the semicircle distribution around zero. 

If the variance of \ogXi is larger than 16, numerical simulation for the central limit is 
difficult. The central limit Y has a big conditional number and all small eigenvalues have 
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Figure 2: The corresponding radius for (p{Y^) and 99(log'^' y) with unity variance of 



big numerical errors. In the further numerical results for cr^ = 64, big eigenvalues of both 
Y and Y^^ larger than unity are numerically very uniformly distributed. 

There are many methods to compare the moments of (j45|) with that of (|24|) to get a sense 
on the similarity between log-semicircle distribution with that of Y. One of the method is 
to find the implied or corresponding radius for each moment by numerically solving 



as a function of k and shown in Fig. [2j The radius are reduced from ri = 2.0816 to 
rioo = 2.0805 for large k with a difference of about 11 x 10~^. The difference of rioo — cq is 
about 3 X 10"^ 

The corresponding radius can also be found by 



that is also shown in Fig. [2l For log Y, the radius rj. are reduced from ri = 2.0817 to 
rioo = 2.0805 for large k with a difference of about 12 x 10"'^. The difference of rioo — co is 
about 6 X 10"^ 

The radius in Fig. [2] shows that log Y can be well-approximated by semicircle distribu- 
tion with radius of i? = Co for cr^ = 1. 

5 Conclusion 

If Xi are positive free random variables with zero mean logXi, central limit for the product 
of xj^^ are bounded between exp(ibco). If logXj has unity variance, cq ~ 2.08 as a 



^ Ji(fcrfc) _ e^/^ 



fcr^. k 



(61) 




(62) 
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special case. If the variance of logXj is far smaller than 16, the central limit distribution is 
very close to log-semicircle distribution. If the variance of logXj is far larger than 16, the 
logarithm of the central limit is very close to uniform distribution. 
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